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We give an overview of recent work on charge degrees of freedom of strongly correlated electrons on 
geometrically frustrated lattices. Special attention is paid to the checkerboard lattice, i.e., the two-dimensional 
version of a pyrochlore lattice and to the kagome lattice. For the checkerboard lattice it is shown that at half 
filling when spin degrees of freedom are neglected and quarter filling when they are included excitations with 
fractional charges ±e/2 may exist. The same holds true for the three-dimensional pyrochlore lattice. In the 
former case the fractional charges are confined. The origin of the weak constant confining force is discussed and 
some similarities to quarks and to string theory are pointed out. For the checkerboard lattice a formulation in 
terms of a compact U(l) gauge theory is described. Furthermore a new kinetic mechanism for ferromagnetism 
at special fillings of a kagome lattice is discussed. 

PACS numbers: 

I. INTRODUCTION 

Charge is quantized in nature. Even subparticles like quarks which can't exist separately (quark confinement) carry a quantized 
charge. It came as a surprise when Su, Schrieffer and Heeger [F2] pointed out that in trans-poly acetylene (CH)„ excitations 
may exist that carry only a rational fraction of the electronic charge provided those polymer chains are properly doped. Since 
then we have got used to the fact that in solids fractionally charged excitations may exist which are either deconfined or confined 
Is'SS]. It has also become clear that two different types of systems with fractionally charged excitations may occur. One type to 
which trans-poly acetylene or graphene [4] belong does not require electron interactions for the existence of fractional charges. 
However lattice degrees of freedom must be included in order that they may form. For the second type of systems electron 
interactions are crucial for the occurrence of fractional charges. The Fractional Quantum Hall Effect (FQHE) fy\ and strongly 
correlated electrons on geometrically frustrated lattices |5] are examples of that category. 

Trans-poly acetylene, graphene and the FQHE are low dimensional systems. Therefore the question arises whether fractionally 
charged excitations may also appear in three dimensions. Indeed, in 2002 it was suggested by one of us 1 5] that excitations with 
charge ±e/2 may exist in pyrochlore and other geometrically frustrated lattices and therefore also in three dimensions. This is 
of interest since in two dimensional systems like the FQH liquid or graphene excitations with fractional charges obey fractional 
statistics [6]. They are neither fermions nor bosons but instead anyons. If fractionally charged excitations would always imply 
fractional or anyonic statistics then one would exclude such excitations in three dimensional systems the reason being that in 3D 
only fermions or bosons can exist as free particles. The possible existence of fractionally charged excitations in a pyrochlore 
lattice excludes a simple one to one correspondence between fractional charges and fractional statistics. 

We argue that a general prerequisite for fractional charges based on electron interactions are strong short-range correlations and 
special lattice fillings. Following the general usage, we will simply call a lattice a "frustrated" one when the short-range elec- 
tronic correlations are incompatible with the lattice structure. It would be more appropriate to speak of "interactions, frustrated 
by the lattice" but that is not customary. 

Here we want to give an overview of recent studies which have been made on c harg e degrees of freedom of strongly correlated 
spinless (or fully spin polarized) fermions on frustrated lattices |5.7. 8 i9lll0llllll2lll3|] . They refer mainly to the pyrochlore lattice 
and its two dimensional projection, i.e., the checkerboard lattice. But also the kagome lattice will be considered in special cases 
|[l4l] . The work on the pyrochlore lattice has been stimulated by experimental findings on LiV204 a transition metal spinel 1 15]. 
It was found that it exhibits heavy quasiparticle behavior, a hallmark of strongly correlated electron systems. Note that spinels 
are of the composition AB2O4 with the B sites forming a pyrochlore lattice. 



II. CHARGES ON FRUSTRATED LATTICES 



In the following we want to concentrate on the crisscrossed checkerboard lattice because a number of features we want to 
point out are more easily visualized on that lattice than on the pyrochlore lattice. The checkerboard lattice can be considered as 
a 2D projection of the 3D pyrochlore lattice (see FIG.IHa), and we treat lattice points connected by a line as nearest neighbors. 
Since we are interested in charge degrees of freedom, we disregard the spin. The following Hamiltonian for fully spin polarized 
electrons or alternatively spinless fermions is assumed to hold 
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FIG. 1: (a) Checkerboard lattice as 2D projection of the 3D pyrochlore lattice fl^. (b), (c) Two examples of allowed configurations on a 
checkerboard lattice at half filling. Occupied sites are connected by thick solid lines as guides to the eye. 




FIG. 2: (a) Adding one particle to the half-filled checkerboard lattice leads to two defects (marked by black triangles) on adjacent crisscrossed 
squares, (b), (c) Two defects with charge ejl can separate by particle hopping without increase in repulsive energy. They are connected by a 
string consisting of an odd number of occupied sites. 



tt = -tY.[c\c^^HA^VY.''^nj . (1) 

{ij) {ij) 

The operators cj create fermions on sites /. The density operators are rii = cjc-. We assume a system of N sites filled with 
N/2 = Y,ni fermions (half filling) and focus on the strong correlation regime, i.e, \t\ <C V. Without loss of generality we assume 

t>0. 



A. Ground-state degeneracy: 

Consider first the case of ^ = 0. The nearest-neighbor repulsions V are minimized if on each crisscrossed square two of the 
sites are occupied while two sites remain empty. This we shall call the tetrahedron rule since it applies equally well to the 
tetrahedra of the pyrochlore structure H7il . All other configurations have a larger potential energy. This implies immediately 

that the ground state is macroscopically degenerate. More precisely the degeneracy is A/deg = (4/3)^^ and therefore the same 
as in the two-dimensional ice model 1 18]. Two of those configurations are shown in FIGs. [Itb,c) where neighboring occupied 
sites have been connected by a solid line. Note that a related Klein type spin model on the checkerboard lattice has been studied 
in Ref. [19.1 and shown to have a similar ground- state degeneracy. It is noticed that each configuration obeying the tetrahedron 
rule, which we call allowed configurations in the following causes a complete loop covering of the plane. When dynamics is 
added to the system the time evolution of the loops will give raise to world sheets in a space-time continuum instead of world 
lines as single particle propagation does. One may speak therefore of a simple form of string theory which is realized here. For 
an extended discussion of that issue see Section IVll 



B. Excitations with charge ±e/2: 

Assume that an electron is added to the otherwise half-filled lattice. In that case two neighboring tetrahedra violate the 
tetrahedron rule because they contain three particles each (see FIG.[2ta). The energy required for adding the particle is 4V. Next 
assume that one particle is hopping to a nearest neighbor site as indicated in FIGs. [2a,b). The new configuration has the same 
total repulsive energy as the old one. This can continue as shown in FIG. Oc) and so on. One notices that the two tetrahedra 
with three particles have separated without increase in repulsive energy. The charge e of the added particle has thus split into 
two fractions e/2. A similar feature is observed if we add an energy V to the ground-state energy. This suffices to break a loop 
as indicated in FIG.[3ta,b). We notice a tetrahedron with three particles and one with one particle only. 



FIG. 3: (a)-(c) Hopping of a particle to a neighboring site: (a) A fractionally charged particle and fractionally charged hole are generated. 
(b)-(c) The two defects (marked by triangles) with charge ±^/2 can separate without creating additional defects and are connected by a string 
consisting of an even number of particles, (d) Example of an allowed configuration on a checkerboard lattice at half filling with possible 
low-order hopping processes [7]. 
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FIG. 4: (a) Spectral function A(k = 7i/2,co) for V = 25t calculated for a \/32 x ^/32 cluster for the effective t-g Hamiltonian ^ with (a) 
g = 0.01 and (b) g = 1. A Lorentzian broadening r{ = 0.1t is used 1^. 

Since a tetrahedron with three particles corresponds to a charge ^/2, the one with one particle only must correspond to a charge 
—e/2. It is an important feature that the fractional charges are always connected by a string of occupied sites, i.e., black lines in 
FIGs. [2land[3l In the case of two charges ^/2,^/2 the string contains always an odd number of sites while for a pair e/2^—e/2 
the number of sites on the connecting string is always even. When dynamics is added to the system the fractional charges will 
separate because each of them will gain kinetic energy. The energy of an added particle with momentum k is therefore 

^(k)=4y + £(ki)+£(k2) (2) 

with k = ki + k2. The form of £(k) is yet unknown. Similarly the energy of a broken string or vacuum fluctuation is 

A^ = y + 8(k)+8(-k) (3) 

where £(k) is the kinetic energy of the fractional charge —e/2. 

Fractionally charged excitations are clearly outside the Landau concept of Fermi liquids. Therefore the spectral function does 
not show a quasiparticle peak like a Fermi liquid does (see FIG.lUa). Instead it shows structures which are partially due to finite 
size effects but have not yet been analyzed in detail. The width of the spectral function is larger than for a particle which is 
prevented from breaking up into two fractional charges e/2. A quasi-particle peak reappears if the ring exchange (see Section 
imi) is so large that confinement (see Section |V]) becomes strong (see FIG. |Hb). A more detailed discussion of the spectral 
function is given in Sect. jV] The remaining question is whether the fractional charges are confined or deconfined. 

C. Side remark: inclusion of spin 

The above considerations referred to spinless fermions. When the spin is included the question arises what happens to it when 
the charge fractionalizes and falls apart into two pieces. To answer it we have to first generalize the Hamiltonian ([TJ to the form 

H = -t'£ (4^;a + H.c.) + y ^ mnj ^uY.ni^riii . (4) 

{iJ)<^ (ij) i 

This is an extended Hubbard Hamiltonian with on-site repulsion U and nt = L^^L^ia- ^^^^ known that in the strong 

o 

correlation limit it can be reduced to a r — / Hamiltonian [|2Q1 of the form 
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H,.j = £ (clc,., + H.C.) + y £ + ^ £ (s,s, - ^) (5) 



where the operators 



C,-a(l -n/-a) 



(6) 



assure that sites are either empty or singly occupied, but never doubly occupied while J = 4t^/U implies an antiferromagnetic 
nearest-neighbor interaction. Consequently the strings which connect the fractional charges represent Heisenberg chains with 
odd (in case of ^/2, ^/2) or even (in case of ^/2, —e/2) number of sites. But a chain with an odd number of sites has a twofold 
degenerate ground state and therefore acts like an effective spin 5" = 1/2. In distinction a chain with an even number of sites has 
a singlet ground state corresponding to 5" = 0. The spin of an added particle is therefore distributed all over the chain and hence 
delocalized over parts of the sample. This a rather unique physical situation. It should be stressed that all of the above features 
are not restricted to the checkerboard lattice but hold for the pyrochlore lattice as well where they are less transparent though. 
A detailed study of the interplay between charge and spin degrees of freedom on the checkerboard lattice at fractional filling 
factors 1/4 and 1/8 can be found in Refs. ||2l||2j23|]. 



III. DYNAMICAL PROCESSES 



The macroscopic degeneracy of the ground state of a half-filled checkerboard- or pyrochlore lattice is lifted when dynamics 
is introduced, i.e., when t with \t\ <C V is taken into account |7]. This lifting takes place to order P /V^, because contribution 
of order t^/V are equal for all configurations. Note that ring hopping processes of order t^/V cancel for fermions. But in order 
P /V^ ring hopping processes connect different configurations with each other. Thus to that order the effective Hamiltonian is of 
the form 



12t^ t t t 

^eff = T- E Ci C ■ C) C ■ C ■ C ■ (7) 

to! 

and the sum is over all hexagons of the lattice. We rewrite this expression in a more pictorial version as 



^eff = -8 L (|<X>)«X>|-|<X>>«^l+H.c.) 

=: -8 L (\b){b\^\B){B\-\a){a\-\a){a\) (8) 

with coupling constant g = I2p /V^ > 0. Again the sum is taken over all, i.e., vertically and horizontally oriented hexagons. The 
symbolic form of //gff is self-explanatory - occupied sites are indicated by a dot, for examples see FIG.Od). The minus signs 
in Eq. ([8]) result from Fermi commutation rules associated with the ring-hopping processes. They depend on the choice of how 
the multi-fermion states on the lattice are enumerated. The simple form of Eq. ([5]) is valid, e.g., for the enumeration used in ?? 
(along the diagonals) In the case of doping //gff is extended to 



Ht-g =H,fi-t'£Pi c]cj + H.c. ) P (9) 

{ij) 

where the projector P projects onto the subspace of configurations with the smallest possible violations of the tetrahedron rule 
compatible with the number of doped particles. For one added particle the number of violations is two. We call the extended 
Hamiltonian the ^ — g model and consider t and g as independent parameters. Thus we shall assume that the two parameters are 
not necessarily restricted to g <c; ^ as required by the strong correlation limit. This is of advantage when numerical calculations 
are performed. The Hamiltonian ([5]) and © does not include 8-site or 10-site hopping processes which become increasingly 
important for increasing ratios t /V . 
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FIG. 5: Height representation for allowed configurations of a \/32 x \/32 checkerboard lattice with periodic boundary conditions at half filling. 
The height field h (numbers in the non-crossed squares) is uniquely defined for a given configuration up to an additive constant M. The field 
f = V/z is indicated by small arrows on the lattice sites. Details of the mapping can be found in the text. The effect on the height fields of a 
ring-exchange process around a hexagon is shown explicitely |7|]. 




(a) (c) 

FIG. 6: (a) Changes in loop topologies by ring hopping processes of Type A (see Eq. (b), (c) Representations of two configurations by 
fully-packed directed loops L11.1 . 

IV. SYMMETRIES AND CONSERVATION LAWS 

As pointed out before, //gff lifts the macroscopic degeneracy of allowed configurations, i.e., those which obey the tetrahedron 
rule. However, the degeneracy is not totally lifted and a residual degeneracy remains. That implies that the total manifold of 
allowed configurations is divided into different subensembles within which all configurations are connected, but which remain 
disconnected by //eff • For their identification it is useful to determine the quantities which remain conserved by //eff. It is easy 
to check that the particle numbers ...,A/4 on the four sublattices into which the checkerboard (or pyrochlore) lattice can be 
divided remain invariant under the operation of //eff. Furthermore, any allowed configuration of the checkerboard lattice can be 
uniquely represented by a vector field f with a vanishing discretized lattice version of curl f. This vector field is obtained for 
the bipartite lattice by assigning an alternating direction of orientation (i.e., clockwise and counter clockwise) to the tetrahedra, 
i.e., crisscrossed squares. Each occupied site is associated with an unit vector in the direction of orientation and each empty site 
with one in opposite direction. This results in a mutual cancellation of the f vectors on a tetrahedron with two occupied and two 
empty sites or more generally when a closed loop is formed. From curl f = it follows that the vector field can be represented 
by f = grad h where the scalar field h defines a hight field up to an arbitrary constant. The height at the upper (right) and lower 
(left) boundary of a finite lattice of A^^ • crisscrossed squares can differ only by an integer — A^y(x) ^ ^y{x) ^ ^y{x) when periodic 
boundary conditions are applied. It is the same for all columns (rows) and defines two topological quantum numbers (Kx,^). 
They are invariant under hexagon hopping processes. As seen from FIG.[5]the application of //gff merely changes the local height 
of two neighboring plain squares by ±2. The doublet (Kj/Nx^Ky/Ny) specifies the global slope of the height field. 

States with (k^, 7^ (0, 0) have a broken lattice symmetry and are charge ordered. For example, when K^/A^x > the charge 
is modulated along a diagonal stripe. Note that violation of the equality Ni = N2 = N3 = N4 also implies a modulation of the 
charge density. 

The different signs in //gff poses a problem when quantum Monte Carlo simulations are applied to the system. Therefore it is 
of considerable help that in important cases this relative sign can be removed by a proper gauge transformation. The system can 
then be transformed to a bosonic one. Without going into details which are found in Ref. ifllll we describe merely the way this 
is achieved. We consider allowed configurations with fixed boundary conditions with an even number of fermions on the four 
boundaries. Open loops are closed as indicated in FIG.Oc) . A given loop covering ^ of the plane is processed as follows. We 
color the background white and alternate the color inside a loop, either from white to dark or reverse. We attach a direction to 
each loop so that the white colored regime is always to the right of the loop (see FIG.[6tb,c)). Next we count the total number 
r of clockwise and / of counter-clockwise loops. The three types of topological changes caused by hexagon flipping processes 
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FIG. 7: Ring hopping changes the number of particles by two on sites marked by blue squares by two. 



with empty center site are shown in FIG.Oa). Note that processes with an occupied center site merely deform loops but do not 
change their topology. The sign change of the processes with an empty center of the hexagon can be cured by transforming the 
configuration to 



1^) . (10) 

The above considerations need not to apply to periodic boundary conditions where only subensembles with even winding 
numbers allow for the coloring prescription given above. Even then it may happen that the coloring is reversed by application of 
^eff while the same loop configuration is recovered. However, from numerical analysis it is found that for tori which preserve the 
bipartiteness of the lattice the transformation applies indeed to the states of lowest energy. In FIG. Oa) we compare the energies 
of a system consisting of 72 sites in which the Fermi sign is taken into account and when it is gauged away, i.e., when the system 
is treated as a bosonic one. The ground-state energy, the first excited states in the (yix^^y) = (0,0) sector and the weights of the 
different configurations are the same in both cases. When (Kx,^) ^ (0,0) we find that in some subensembles the ground- state 
energies are higher for fermions than they are for bosons, i.e., when the fermionic sign is removed by a gauge transformation. 
All the considerations are limited though to hexagon ring hopping processes and do not apply to higher-order processes. 

It should be also mentioned that the sign of g is irrelevant. By multiplying all configurations by a factor where is the 
number of fermions on the sublattice in FIG. H the sign of all matrix elements of ^eff is changed. Invariance under this gauge 
transformation demonstrates a global —g symmetry of the spectrum. 



V. CONFINEMENT VS. DECONFINEMENT 



The question of confinement of fractional charges is closely related to the form of the ground state of the half-filled lattice. 
For example, when the ground state is charge ordered then one expects the fractional charges of a doped system to be confined, 
because a separation of them results in disorder of the ground state. Therefore we discuss first the character of the ground-state 
wavefunction. The Hamiltonian is given by Eq. ([5]) but we want to supplement it following Rokhsar and Kivelson [24] by an 
extra term which counts the number of flippable hexagons. Depending on its sign, configurations with such hexagons are either 
favoured or partially suppressed. The Hamiltonian which we consider reads therefore 

= /^eff+^i:{^,Q}(|<X>)(<%>l + l<X>)«^l) (11) 
where the sum is over all flippable hexagons with occupied or empty center. Two limiting cases are particularly simple: 

(i) [A +oo: all configurations without flippable hexagons (frozen configurations) are ground states with 
energy £^ = 0. 

(ii) jLi^ configurations with maximal number of flippable hexagons Nfi become ground states. Numerical 
Monte Carlo calculations for checkerboard lattices up to 1000 sites yield configurations with a "squiggle" 
type structure 1 10] or slight variations of it as ground states. The results suggest that in the thermodynamic 
limit they are all in the (K;^,^^;) = (0,0) subspace. Due to the large unit cell of the squiggle configuration 
and the conservation laws the number of different ground states is ten. 
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FIG. 8: (a) Ground- state energy and energies of the lowest excited states of the effective Hamiltonian //gff in subspaces with different global 
slopes (Kjc, Ky) = (Kjc, 0) for a 72-site cluster. Right side: Same system, but assuming same signs for all matrix elements ("bosonic calculation"), 
(b) Energies of the ground state and lowest excited states for all subensembles of a 72-site half-filled checkerboard cluster for different values 
of /i of the g-}A Hamiltonian. Level crossing of ground states occurs only at ^ = g. The insets indicate different phases: Maximal flippable 
plus fluctuations for /i < g, a critical point id = g where the ground state is an equally weighted superposition of all configurations, and frozen 
configurations as ground states for /n> 1.11.1 . 

(iii) /i = g > 0: This is a particularly interesting case because it is exactly solvable ll24l] . The ground states have 
energy E = and some of them are liquid like. By using the gauge transformation ([TOl) we change the sign 
of the second term in ^ and rewrite Hg=^ in the form 

Hs=, = g I [(l<X>)-|<^))x(«%>|-(<X>l)] (12) 

which is a sum over projectors. Therefore all eigenvalues must be non-negative. For each subensemble / 
the ground-state wavefunction is then given by an equally weighted superposition of all connected config- 
urations 1 4^^), i.e., 

|v)/«)=A£|c«) (13) 

i 

where A is a normalization pref actor. Using the Perron and Frobenius Theorem (see below), it follows 
that these are the unique ground states in each subensemble. The coherent superpositions are resonating 
Valence Bond (RVB) states of the form originally discussed in |25.26]. We want to draw attention that 
a RVB state is obtained only when a gauge transformation of the form ([TOl) exists which ensures that all 
off-diagonal matrix elements have the same sign. Otherwise the energy is most likely larger than zero. 

Next we explore the phase diagram as function of ^/g by exact diagonalization of Hgjj for clusters up to 72 sites. The details 
are found in Ref. ifToll so that we state here merely the results: 

(i) For in> gwe find the same features as for /i ^ oo, i.e., all frozen configurations are ground states. 
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FIG. 9: Half filling: (a) Charge density distribution for one of the two ground states, (b) Corresponding density-density correlation function 



absolute value. Red or blue color represents a positive and negative value, respectively, (c) Classical density-density correlation function. 



(ii) For /i < g we find two ground states in subensembles with (k^, k^) = (0,0). They are (M, ...,A^4) = (6, 6, 12, 
12) and (12, 12, 6, 6) respectively. They consist of superposition of configurations with maximal number of 
flippable hexagons and of fluctuations around these configurations. We expect that in the thermodynamic 
limit the 10-fold degenerate squiggle phase instead of a two-fold degenerate ground state is recovered. The 
average weight of configurations with maximal number of flippable hexagons decreases as ju increases until 
at /i = g all configurations of a subensemble have the same weight. 

The energy of the ground state and of the lowest excited states as function of /i/g is shown in FIG. [8] for a cluster of 72 sites. 

It is interesting to calculate the density-density correlation function for a 7 2- site cluster at /i = when one of the two degenerate 
ground states is used for \|/o, i.e.. 



Qo = (¥o I nm I \|/o) - (\|/o I m I \|/o)(\i/o I m \ \|/o) 



(14) 



Here the subscript denotes a site near the center of the cluster. The charge order for the state M , A/4 = (6, 6, 12, 12) yields 
stripes of average occupation 1/3 and 2/3 (see FIG.[9ta). The result for Qo is shown in FIG.[9tb). It differs only slightly from 
the corresponding classical correlation function plotted in FIG. [^c) which is obtained by setting ^ = and summing over all 
degenerate configurations. 

We are now in the position to discuss the problem of confinement of fractional charges. It is obvious that for in> g, i.e., 
when the ground state energy is £" = it does not matter how far two charges e/2 of an added particle are separated from each 
other since the energy is always 4V. However that is different for ^ = 0, the case of physical interest. The energy change can be 
decomposed into local contributions which include all hexagon hopping processes involving site /. As seen in FIG. fTOla-c) 
Ei increases along the string connection between the two fractional charges (compare with FIG.O. This is plausible since ring 
hopping requires alternating empty and occupied sites. In the vicinity of a string of occupied sites these processes are reduced 
and increases. The same holds true when a pair e/2, —6/2 is separated. Therefore in both cases the fractional charges are 
confined by a constant confining force, similarly as quarks I27ll28|] . The energy increase is linear with charge separation. From 
the numerical results we can deduce an energy increase of AE ^ 0.2\ g\r where r is the distance of the fractional charges in units 
of the site distance a. The energy increase corresponds to a string tension T = 0.2\g\.lt should be noticed that the restoring force 
is weak. Assume that V = lOt in which case g = 0.12 and AE = 0.024ta. The two fractionally charged particles e/2 will form 
a bound state with a radius of order 50 - 100 lattice distances. Thus as soon as there is a small doping concentration the average 
distance between the fractional charges will be smaller than the diameter of the pairs and a plasma will form. Note that when r 
exceeds a critical value so that AE = 0.2g • Vc then it is favourable to create fractionally charged particle-hole pairs e/2, —e/2 
in order to reduce a further energy increase with r. Again, this resembles the generations of quark- antiquark pairs in form of 
pions, when a quark is separated too far from the other ones. 

One would like to know how fractionalization of charge could show up in experiments. One quantity which depends strongly 
on it is the spectral function A (k, co) = (k, co) + A~ (k, co) which is defined through 



A+(k,co) = lim --Im(\|/Q^ 

r|^0+ n \ 



1 



CO + /r| + £0 - ^ 



A (k,co) = lim --Im(\|/^ 
r|^0+ 71 \ 



1 



co + /r| -^0 + ^ 



(15) 



Here |v)/^) is the ground state with energy £0 of the half-filled lattice. We have calculated A (k,(o) by using the effective Ht-g 
Hamiltonian for a finite cluster and different values of g (see FIG.|4li. For small values of g the diameter of the bound pair (e/2, 
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FIG. 10: (a)-(c) Local loss of kinetic energy due to the separation of two (static) fractionally charged defects, i.e., particles or holes (fcp's 
marked by light red squares or fch's marked by dark blue squares). The radii of the circles are proportional to the local energy loss, (d)-(f) 
Red (blue) circles show an increase (decrease) of the local density (vacuum polarization due to the two fcp's or the fcp-fch pair). 



e/2) formed by an added particle exceeds the cluster size. Therefore the fractional charges appear to be deconfined and as seen 
in FIG. Ufa) there is no quasiparticle peak in A(k,co). This is different when g is sufficiently large, e.g., for g = 1. Here the 
bound pair remains within the finite cluster and therefore a quasiparticle peak does appear (see FIG.Ill^b). The internal degrees 
of freedom of the bound (e/2, e/2) pair yield additional structure in A(k, co). Note that in FIG. Ufa) as well asH^b) the bandwidth 
of A(k, co) is nearly twice as large than as without charge fractionalization in which case it would by 8t. 

Hand in hand with the separation of two fractional charges goes a polarization of the vacuum between the two objects. The 
modifications of the vacuum in the vicinity of the string connecting the two fractional charges is shown in FIG.fTOld-f). We may 
consider the vacuum polarization as origin of the confining force. 

The above considerations apply to the checkerboard lattice. As far as the pyrochlore lattice is concerned we do not yet 
know whether fractionally charged particles are confined or deconfined. There are numerical indications that the charges are 
deconfined here due to the smaller effect of a string of occupied sites on ring hopping. But only detailed numerical work can 
provide a solid answer [29^] . 



VI. STRINGS AND THEIR TIME EVOLUTION 





FIG. 1 1 : Continuous representation of loop dynamics due to Hgff . (a) Time evolution of loo^s due to B processes which conserve the topology, 
(b) - (d) The same for A processes where Hgff induces three kinds of topological changes fl3|]. 

As it was shown above, when spinless fermions occupy a checkerboard lattice at half filling the strong correlations result in 
a complete loop covering of the plane. The time evolution of these loops results in world sheets which replace world lines of 
single particle propagation. Therefore we are dealing here with a simple form of a string theory, a very active and challenging 
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topic of present days field theory lll2ll3Q|] . It is interesting to consider the loop dynamics due to ^eff in a continuum limit. This 
is shown in FIG.[TT] There are processes denoted by B in Section Hill which change only the shapes of the loops continuously. 
They are visualized in FIG.[Tna). In distinction A processes in ifgff change the topology of a given loop covering. They were 
shown in FIG.Oa) and give raise to a time evolution shown in FIG. fTTl b-d). 




FIG. 12: When an energy A£ > V is added to the system a loop is broken and an open string is generated. Thereby one end of the string is 
touching a closed loop, (a) Time evolution in a continuous representation, (b) Time evolution of two (^/2, —e/2) pairs and formation of e and 
—e particles. 

When an energy AE" > V is added to the ground state a loop can break up and a pair of charges — is generated. This 
process is shown in FIG.[T2ta). When two or more pairs are generated new combinations e/2, e/2 and —e/2, —e/2 may form. 
They correspond to the creation of new particles out of the vacuum (ground state). Note that in order that, e.g., the two fractional 
charges e/2, e/2 combine to a particle with charge e they must be situated on different sublattices. This is schematically shown 
in FIG.fT2tb). Fractionally charged particles are always situated at the end of open strings. In case of ^/2 (but not —e/2) the 
open end is always touching a closed loop. Similar considerations hold for pyrochlore lattices at half fillings. 

VII. DERIVATION OF A GAUGE THEORY 

The problem of confinement vs. deconfinement of fractional charges on a checkerboard lattice was described in Sect. El It 
can also be discussed in terms of a field theory, more specifically a gauge theory. Since the effective Hamiltonian conserves the 
number of particles on a crisscrossed square the theory is invariant with respect to gauge changes on those squares lfT2|] . This 
enables us to rewrite the Hamiltonian in form of a U(l) lattice gauge theory and to compare it with that of other models. Similar 
theories were derived before for the quantum-dimer model [3132,1 , for a three dimensional spin system [ 3334] and for the ice 
model 1I35I1 . The derivation follows the one given by Polyakov ll36ll37ll for compact quantum electrodynamics. By connecting 
the centers of crisscrossed squares we obtain a square lattice with the particles sitting on links instead of sites. Define for each 
link X, X + between neighboring lattice sites x and x + e^- (j = 1,2) a variable /tj(x) with integer eigenvalues. In order to express 
the effective Hamiltonian in terms of /tj(x) we first introduce its canonical conjugate, i.e., the phase ^j(x) G [0,7c]. We note that 
exp. [±/^j(x)] acts like a ladder operator. In rewriting ring hopping in terms of /2j(x) and $j(x) we must ensure that a link is 
either occupied by one particle or unoccupied. That leads to the following form of the effective Hamiltonian ^ 




where the first term limits the eigenvalues of ^j(x) to and 1 while the second term describes ring hopping. Note that here 
hexagon ring hopping goes over into hopping around two neighboring squares (double plaquettes) and that the signs of the 
phases alternate around the polygons. We have assumed that by means of ([TOl) the sign change in ^ has been removed. 
Next we introduce staggered gauge and electric fields on the biparticle lattice 

A,-(x) = (-ir'+^^$;(x) 

£,-(x) = . (17) 

The theory must incorporate the constraint that each lattice site is touched by exactly two occupied links. This follows from the 
rule that each tetrahedron contains two particles. Here this constraint reads 



(A,£, (x)-p(x)) |Phys)=0 



(18) 
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with the lattice divergence defined by 

AjEj{x)=Ei{x)-Ei{x-ei)^E2{x)-E2{x-e2) . (19) 

It is noticed that the constraint has here the form of Gauss' law. In terms of Ej{x) and Aj{x) the effective Hamiltonian ([T6b 
becomes 

//eff = liint/EfeW-7)-2gIcos(£Aax)+£Af(x-e,)) . (20) 

°° xj \ ^ / xj \ n □ / 

The oriented sum of the vector potential around one plaquette is 

^A^(x) = Ai(x) -Ai(x + e2) -A2(x) +A2(x + ei). (21) 

□ 

The Hamiltonian resembles the one of compact quantum electrodynamics in 2+1 dimensions ll36ll37ll but is not identical with 
it. It is interesting that in the latter case two charges are confined and that the energy grows linearly with there distance as it 
is the case in our model. However, despite of the above mentioned resemblance the physical origin of confinement is different 
in the two cases. As it turns out the U(l) gauge theory formulation of our model system is gratifying but does not bring much 
additional insight. 



VIII. FERROMAGNETISM GENERATED BY KINETIC PROCESSES 



The most common origin of ferromagnetism is spin exchange between electrons. The latter may belong either to different 
atomic sites or to the same site where intra-atomic exchange is the origin of Hund's rules. Pauli's principle forbids electrons with 
parallel spins to come too close to each other and reduces this way the mutual Coulomb repulsion. This should be compared with 
superexchange, the standard mechanism for antiferromagnetism. Here it is the kinetic energy which is optimized in the magnetic 
phase. Hence it is in general a competition between an optimization of the Coulomb repulsion and of the kinetic energy which 
favours magnetic order. In passing we mention other sources of ferromagnetism such as RKKY interactions in metals or double 
exchange, to name of few. 




FIG. 13: Panels (a) and (b) show two different allowed configurations on a kagome lattice which fulfill the constraint of zero or one electron 
per site and one electron of arbitrary spin per triangle. The arrows indicate possible ring-hopping processes. All allowed configurations can be 
represented by colored dimer coverings on a honeycomb lattice shown in (c) and (d) which is obtained by connecting the centers of triangles. 
A blue dimer correspond to a spin up and a red dimer to a spin down particle. 

It might appear surprising that ferromagnetism may also be caused by purely kinetic processes. We know of two specific 
examples where this is the case. One is the ferromagnetic ground state discovered by Nagaoka | 38]. It is due to a single hole 
moving in an otherwise half-filled Hubbard system in the limit of infinite on-site repulsion U. The proof is based in a theorem 
due to Perron and Frobenius. That theorem states (see, e.g., Ref. [39] ) the largest eigenvalue of a symmetric n • n matrix with 
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only positive matrix elements is positive and non-degenerate, while the corresponding eigenvector is nodeless, i.e., we may 
choose it to have only positive components. The theorem applies only to systems with a finite-dimensional Hilbert space. The 
same theorem is the basis of three-particle ring exchange in ^He. In both cases ferromagnetism results from the motion of 
fermions because the ground state wavefunction is the smoothest in this case (it is nodeless) and has the lowest kinetic energy. 
Here we want to point out another source of ferromagnetism based on purely kinetic effects. We demonstrate it by considering 
a partially filled kagome lattice with electrons described by an extended Hubbard Hamiltonian (|4]). We focus on the case of 1/6 
filling implying one electron per triangle. In the limit of ^ oo double occupancies of sites are excluded and strong correlations, 
i.e., 1^1 <C y are assumed. When t = the ground state is macroscopically degenerate since all configurations with one electron 
of arbitrary spin orientation per triangle are ground states. For examples see FIG. [T3ta,b). By connecting the centers of the 
triangles of the kagome lattice we obtain as the medial lattice a honeycomb lattice. Here electrons sit on links instead of lattice 
sites. The different ground-state configurations correspond here to two-colored (spin) dimer configurations (see FIG. [T3lc,d). 
They are orthogonal because possible wavefunction overlaps are neglected. 




-0.2 0.2 0.4 0.6 -0.2 0.2 0.4 0.6 

J/g J/g 



FIG. 14: Exact diagonalization of the two-color dimer model on 24- site honeycomb cluster. The upper panels show the ground- state energies 
of different sectors as a function of next-nearest neighbor coupling J/g. The lower ones show the expectation values of the spin part of the 
Hamiltonian. The ground state is denoted by |0) lH^ . 

When we allow for ^ 7^ the ground-state degeneracy is lifted. To lowest order of degenerate perturbations theory the effective 
Hamiltonian acting on the subspace of configurations with precisely one dimer touching each site is written as 

//hex = -8 L (lOXOl + lOXOl +H.C.) (22) 
{oHaB.} 

where g = 6p/V'^. The sum is overall hexagon and spin combinations symbolized by the three marks. Note that ^hex causes no 
fermionic sign problem. If we would ignore the spin the model would go over into the quantum dimer model (QDM) studied in 
Ref . ||4Q|i . The ground state of the QDM is known to be three-fold degenerate in the thermodynamic limit. It corresponds to the 
valence-bond solid (VBS) plaquette phase with broken translational symmetry. Here we are interested in spin correlations and 
notice that //hex conserves both the total spin S as well as its component S^. Without proof which is found in Ref. [14] we state 
that the Perron-Frobenius Theorem is applicable to the present system. 

Returning to //hex we remark that by a simple gauge transformation the sign of the plaquette flip in Eq. ([22l) can be always 
chosen negative, irrespective of the sign of t 1 14]. We therefore choose all off-diagonal matrix elements of //hex to be non- 
positive. But before applying the Perron-Frobenius Theorem we have to discuss the problem of ergodicity. The Hilbert space 
under consideration is broken into different sectors corresponding to different which are not connected by //hex- Therefore 
each sector must be considered separately. For the S^^^ = 5'max = Ne/2 sector, where Ne is the number of electrons the ground 
state is unique and fully spin polarized. The situation is more complex for the other sectors with S^^^ 7^ 5'max. Although it can 
be shown that a fully spin polarized state is also a ground state in those sectors, it has not been possible to prove convincingly 
ergodicity, i.e., that the fully spin-polarized ground state is not only one among others. However, the numerical studies on 
clusters strongly suggest that the ground state is indeed unique. 

In order to study the robustness of the kinetic ferromagnetism, we introduce an additional next-nearest neighbor antiferromag- 
netic interaction so that the total effective Hamiltonian becomes //eff = //hex + //spin with 
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Hspm=J X (SiSj-jmrijj . (23) 

Numerical diagonalization of clusters up to 24 sites show that the ground state remains fully spin polarized up to values of 
J/g < {J/g)c ^ 0-2. This is shown in FIG. fT4l which demonstrates a considerable robustness of kinetic ferromagnetism. 

A ferromagnetic ground state is also found for a filling factor of 1/3. The origin of kinetic ferromagnetism in the present 
examples differs from the ring exchange proposed by Thouless for ^He [41] insofar as there the particles remain at their orig- 
inal location when they cyclically permute while here they actually move around. It differs also from the so-called flat band 
ferromagnetism of Mielke which he predicted for a kagome lattice with fillings between 5/6 and 1 1/12 ll42l43h . In his case any 
value ofU > is sufficient and V = while in our case U ^ oo and V/ms large. In Mielke's case a completely flat band in 
a single particle description of electrons on a kagome lattice influences so strongly the Stoner criterion for the occurrence of 
ferromagnetism that any value ofU >Q leads to a ferromagnetic instability. 
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